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ABSTRACT 



The efficient use of a search-surface radar or sonar, in 
which one or more targets appear on the screen intermittently 
usually demands a device for tracking the targets auto- 
matically. Such a device, called a "track while scan system", 
must make an estimate of each target's instantaneous position 
from the sampled-data information provided by the radar. 

For this purpose, an ct-6 filter and an optimal Kalman 
filter, that must track maneuvering targets, are analyzed 
here and com.pared in terms of tracking accuracy for tactical 
applications . 
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I. 



INTRODUCTION 



Systems such as search-surface radar and sonar, in which 
the input data arrives intermittently, frequently require a 
device for continuously estimating the "present" value of 
the input. In radar terminology, this device is called a 
"track while scan system" . 

In a more general sense, the term "track while scan" 
system, m.ay denote any system, which estimates the "present" 
value of a signal from the "past" sampled values of the signal 
the sampling taking place at regular intervals. 

Vehicles, without maneuvering, of the class under consider 
ation (such as aircraft, ships, and submarines) generally 
follow straight line constant velocity trajectories. If the 
vehicles were not able to deviate from these trajectories, 
i.e., could not maneuver, then the tracking problem could be 
solved quickly and sim.ply using standard filtering such as 
an a-B filter. 

Historically, a-$ filters were designed to minimize the 
mean square error in filtered position and velocity under the 
assumption of small velocity changes between data samples. 
Thus, most a-3 filters have little capacity to track targets 
that either accelerate or maneuver (changing direction) . An 
important early paper [1] defines 6 in terms of a with a a 
design parameter. This greatly simplified the use of maneuver 
detectors by reducing the design problem to a single variable. 
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A later advance in a-B filter design was achieved by 
minimizing the mean square error in predicted position 
[ 2 ] . 

A system that developed using the above technique is the 
VEGA LN I (Thomson-CSF, France) . 

The final choice of parameters for the a-B tracker will 
always end in a compromise between smoothing the input noise 
(measurement errors) and retaining some ability to follow a 
maneuvering target. Various performance measures have been 
used for this compromise. For example, steady state noise 
variance reduction and the ability to follow a target capable 
of responding to impulse accelerations are used as perform.ance 
measures in [1] to derive good steady state filter parameters. 
Time varying noise variance reduction and the ability to follow 
a randomly maneuvering target are performance criteria which 
follow from [3] . In this paper, such a filter is analyzed, 
the frequency response, stability, noise characteristics and 
transient error are derived and plotted, and a scheme for 
optimizing the two dynamic parameters (a-B) is suggested. 

A criterion for a-B filter design with a numerical and 
graphical example is also given. 

Since the majority of tactical weapons systems requires 
that manned maneuverable vehicles, such as aircraft, ships, 
and submarines, be tracked accurately, an optim.al Kalman filter 
has been derived for this purpose, based on the early work of 
[4,5]. The target model for tracking applications must be 
sufficiently simple to permit ready implementation in weapons 
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systems for which computation time is at a premium yet 
sufficiently sophisticated to provide satisfactory tracking 
accuracy . 

The target acceleration selected for the Kalman filter, 
and hence the target m.aneuver is correlated in time; namely, 
if a target is accelerating at time t, it is likely to be 
accelerating at time t+x for sufficiently small x. By sim- 
plifying the maneuvering model used in the Kalman filter above, 
the state vector can be reduced from six to four elements. 

The model simplification - simplified Kalman filter - is 
achieved by assuming (incorrectly) that the vehicle's change 
in velocity is uncorrelated between samples; i.e., the maneuver 
is white. In the comparison which follows, the Kalman filter 
and Sim.plified Kalm.an filter are compared with the a-g filter 
in a variety of tactical environments with tracking sensors, 
utilizing Monte Carlo simulation techniques on realistic 
target trajectories to verify their theoretical performance. 
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II . FUNDAMENTALS OF TRACK -WHILE -SCAN 

Any system which performs a tracking function must obtain 
and utilize the basic target parameters of position and rates 
of motion. In the earlier system, both position data and 
velocity data were used to maintain the tracking antenna on 
the target at all times, thus limiting the system to the classi- 
cal one target-at-a-time tracking function. 

In a track-while-scan system, target position must be 
extracted and velocities calculated for many targets without 
holding the radar antenna fixed on one target. Obviously in a 
system of this type, target data is not continuously available 
for each target track at a rate dependent upon the scan rate 
of the system. 

In a typical TWS system the data rate is one unit of data 
per second or a scan rate for the search radar antenna. 

Since the antenna is continuing to scan, some means of 
storing and analyzing target data from one update to the next 
and beyond is necessary. The digital computer with its memory 
and computational capability is employed to perform this 
function and also: 

(i) To provide a tactical picture 
A display of the tracks of all vehicles observed by sensors 
showing present positions, courses and speeds etc., is essen- 
tial for the deployment of a ship and its weapons. The computer 
enables rapid use to be made of sensor information thereby 
ensuring that the picture is accurate and up-to-date. 



(ii) For use in processes such as threat evaluation 
and weapon assignment 

The computer can forecast future likely positions of tracked 
vehicles and rapidly perform necessary calculations to assist 
operators in the assessment of threats and in the optimum 
weapon deployment to deal with them. 

(iii) To provide target information for weapon 
deployment 

The computer can be used to generate smooth tracks from noisy 
information thereby being able to pass information to a weapon 
sensor more accurately than information derived from a single 
plot on display. Hence, the central concept underlying any 
TWS system is that the sensor itself continues to perform its 
primary function of search (scanning) and data input while the 
remainder of the system performs the target tracking function. 
The sensor function simply provides target position data to 
the com.puter subsystem where target velocities and position 
prediction are calculated. In a military application the 
major advantage of a TWS system is the elimination of the 
process of target designation from a search radar to a fire 
control radar. 

The tracking information, developed in the TWS system, 
is used as a direct data input to the computation of a fire 
control solution. 

Therefore, as soon as a target is detected a fire control 
solution is available without the inherent delay caused by 
the designation process. The time required from first detection 
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to fire control solution is on the order of inilli-seconds for 
a TWS system as opposed to seconds or even minutes for a 
manually designated system employing separate search and fire 
control sensors. 

The focus of the following chapter has been to answer the 
question: "What functions should the TWS system perform in 

order to com.bine the search and tracking tasks into one 
integrated unit?" 

A. A TRACK WHILE SCAN METHOD 

The method of solving the track while scan problem is 
based upon the assumption that the radar furnishes target 
position information once each scan. The scheme can be imple- 
mented by a combination of special radar circuits (hardware) 
and software. Existing systems also provide for operator 
modification of system tracking. 

Any track-while-scan system must provide for each of the 
following functions: 

(i) Target detection 

(ii) Generation of target acquisition and tracking 
"windows" or "gates" 

(iii) Target track initiation (assignemnt of targets 
to track files) 

(iv) Target data input and track correlation 

(v) Track "window" prediction, smoothing and 
positioning 
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1. Target Detection 



Target detection, localization and designation are 
accomplished by the radar sub-system in the usual sense. 

Hence, the TWS receives the following data from the outside, 
via wire links: 

(i) The surveillance antenna azimuth 

(ii) The radar synchronization (presyne) 

(iii) The surveillance radar video 
2 . Generation of Target Acquisition and Tracking "Windows” 

A "window" can be defined as a sm.all volume of space 
initially centered on a target, which will be monitored on 
each scan for the presence of target information. Each initial 
target detection will cause a relatively large acquisition 
"window" to be generated, centered on the position of the 
detection. When a target is initially detected the algorithm 
receives only the position data for that initial, instantaneous 
target position. The acquisition window is then generated, 
for example: 

Initial position = (Range, Bearing, Elevation) 

Range window = R ± 1000 yars (914.1 meters) 

Bearing window = B ± 5° (0.0873 radians) 

Elevation window = E ± 5° (0.0873 radians) 

The acquisition window. Fig. la, is large in order to allow 
for target displacement during the following scans of the 
radar. If on the next radar scan the target is within the 
acquisition window, a tracking window is generated in the 
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Azirnuth Window 
10 degrees 
0.1745 radians 





b. TRACKING VTIITO-J 



FIG. 1. 



TRACK WHILE SCAN VOLUMETRIC WINDOWS 



I 




same manner as the acquisition window. Although Fig. lb 
shows only the very small, 120 yards, 1.5°, 1.5° tracking 
window, in actual practice intermediate window sizes are 
generated until a smooth track is achieved. 

3 . Track Initiation 

Concurrent with the generation of the acquisition 
window a track file is generated in order to store the posi- 
tion and window data for each track. In addition to the basic 
position and window data, calculated target velocities and 
accelerations are also stored in each track file. Track files 
are stored within the digital computer (or processor) sub- 
systems memory and the data is used to perform, the various 
calculations necessary to maintain the track. 

Each track file occupies a discrete position of the 
digital computer's (or processor's) high speed memory. As 
data are needed for computation or new data are to be stored, 
the portion of memory which is allocated for the required data 
will be accessed by the system programs (software) . In this 
manner a diversity of data in addition to the tracking data 
may be stored in the "track" file, for example, ESM data, IFF 
inform.ation • The generation of the track file begins with the 
initial storage of position data along with a code to indicate 
that an acquisition window has been established. If target 
position data is obtained on subsequent scans of the radar, the 
file is updated with the coordinates; the velocities and 
accelerations are computed and stored, and the acquisition 
window code is canceled. 
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The acquisition window is then decreased in size to 
the tracking window and the track code is stored which indi- 
cates an active track file. 

As the radar continues to scan, each input of data 
is compared with the window positions of active track files 
until the proper file is found and updated. However, it should 
be noted that the search for the proper track file is generally 
not a sequential one-to-one comparison. This method is much 
too slow to be used in a system where speed of operation is 
one of the primary goals. 

This idea of comparing output data with window posi- 
tions leads us to the problem of correlation of data input to 
track files and what to do if correlation is ambiguous. 

4 . Resolution of Track Ambiguity 

Track ambiguity arises when either multiple targets 
appear within a single track window or tv/o or more windows 
overlap on a single target. 

This occurrence can cause the system to generate 
erroneous tracking data and ultimately lose the ability to 
maintain a meaningful track. If the system is designed so 
that an operator initiates the track and monitors its progress, 
the solution is simply for the operator to cancel the erroneous 
track and initiate a new one. 

For systems which are to be automatic, software (pro- 
gramming) decision rules must be established which will enable 
the tracking program to maintain accurate track files. Decision 
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rules as to target input data and track correlation can be 
logically developed with an understanding of the definition 
of the track ambiguity problem. 

5 . Track Window Prediction, Smoothing, and Positioning 

In a track-while-scan system tracking errors also 
exist due to target motion. The tracking window now has 
replaced the "tracking antenna" and this window must be posi- 
tioned dynamically on the target in a similar manner as was 
the "tracking antenna". However, there is no "servo" system 
to reposition and sm.ooth the tracking window's motion. This 
repositioning and smoothing m.ust be done mathematically within 
the TWS algorithm. To this end, smoothing and prediction equa- 
tions (Eqs. 2. 1-2. 3) are employed to calculate the changing 
position of the tracking window. Instead of the system "lagging" 
the target the tracking window is made to "lead" the target 
and smoothing is accomplished by comparing predicted parameters 
with observed (measured) parameters and making adjustments 
based upon the errors derived from this comparison. 

Figs. 2 a-b and 3, illustrate the TWS principle, 
track positioning and the TWS general processing loop of the 
alpha-beta type, respectively. 

B. DEFINITION OF TRACKING EQUATIONS FOR a~B FILTER 

Tracking in a track-while-scan (TWS) system consists 
basically of two functions: 

a. Smoothing 

b. Correlation 
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FIG. 3. TWS GENERAL PROCESSING 



Smoothing is the processing of the sensor reports to derive 
an estimate of both target velocity and position. 

Correlation is the sorting of sensor reports into groups, 
to determine which belong to which target then predicting a 
new coordinate on which to center the correlation region. 

The constants of proportionality, a and 3/ used in 
correcting the position and velocity, respectively, of the 
estimated target course, completely characterize the perform- 
ance of the TWS system. These constants are the dynamic 
parameters or so-called "sm.oothing constants" of the system. 

Finally, the simplest case target tracks are based on 
smoothing and prediction of an alpha (a) - beta ( 3 ) tracker 
operating in a cartesian coordinate reference frame. The 
information is ordinarily obtained from a coordinate converter 
operating on the raw polar to cartesian transform. The a-3 
filter described by the well-known a-3 tracking equations [6,7] 

Xg = Xp + a{X^ - Xp) smoothing Equation (2.1) 

3 (X^ — X^) 

^ prediction Equation (2.2) 

X^ = X^~^ + T prediction Equation (2.3) 

IT ^ O 

where : 

= smoothed position 
Xp = predicted position 
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= smoothed velocity 
measured position 

filter param.eters or filter gains 
sampling time or time between detections 

Substituting Eq. (2.3) into Eqs . (2.1) and (2.2) yields 



M 
a , 6 

T 



N N-1 N-1 N 

Xg = (l-a)Xg + (1-a) TVg + axj^ 



(2.4) 



— — — X^~ ^ + ( 1 — ft ^ V^~ ^ + — X^ 
Vg - T Xg + (1 6)Vg + T 



(2.5) 



x: 



,N+1 



V- T 



( 2 . 6 ) 



or 



and 



^s’ 


N 


(1-a) 


(l-a)T 




"7 






6 

T 


(1-6) 










N 



T 



N-1 



(2.7) 



N+1 



N 



[Xp] 



= [1 T] 



X, 



( 2 . 8 ) 



V, 
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These equations represented in block diagram form appear in 
Fig. 4. 




Fig. 4. Block Diagram of an a-3 Tracker 



where the following matrices and vectors are defined: 





(1-a) (l-a)T 




a 


A = 




B = 






-6/T (1-g) 




B/T 



C = [1 T] X = 



1 . Stability Analysis 

Since the system characteristic equations are in the 
form of difference equations, Z-transform theory is helpful 
in determining filter transfer functions and stability criteria. 
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For the system of the form; 



- c='n> 

or 

Vl ' <* - ^ 



thus , 



G(Z) 



X(Z) 

X„(Z) 



= [ZI - (A-BC) ] B 



(2.9) 



Using Eq. (2.9), the Z-transform for each of the several 
components of Eg. (2.7) and (2.8) has been computed and is 
given below. Hence, the a-B tracker transfer functions; 



Ss^2) 



G 



VS 



(Z) 



G^p(Z) 



G 



VP 



(Z) 



_ Z [g (Z-1) + 3] 



X„(Z) 


D(Z) 


Vg(Z) 


|z(Z-l) 


X„(Z) 


D(Z) 


Xp(Z) 


g (Z-1) + 6Z 


Xm(Z) 


D(Z) 


Vp(Z) 


|(z-l) 


X„(Z) 


D(Z) 



Smoothed (2.10) 

position 



Smoothed (2.11) 

velocity 



Predicted (2.12) 
position 



Predicted (2.13) 
velocity 
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where the characteristic equation of the system is 



D(Z) = I ZI - a| = 0 



or 



D(Z) = Z^ - (2-a-6)Z + (1 - a) = 0 (2.14) 

The stability of the system can be determined from the loca- 
tion of the roots of the characteristic equation above. 

A few methods are available for determining whether 
or not a polynomial in Z contains a root or roots on or out- 
side the unit circle. One method is to modify the Routh- 
Hurwitz stability criterion. The Routh stability criterion 
tells whether or not any of the roots of a polynomial lie in 
the right half of the complex plane. 

Since the following transformation 



maps the interior of the unit circle in the Z plane to the 
lef-half r plane, with this transformation, permits the 
application of the Routh stability criterion to the polynomial 
in r as in continous-time systems. 

Now Eq. (2.14) becomes a polynomial in r, in the 

form 
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2a 



6 



0 



(2.15) 



6r + 2ar + 4 - 



So, the resulting requirements for stability are 



a > 0 , 6 > 0 and (2a + B) < 4 (2.16) 



An additional stable condition exists when 6=0. Thus the 
resulting necessary and sufficient conditions for the stability 
of the track-while-scan system are 

a > 0 , 6^0, and (2a+6) £4 (2.17) 

These inequalities determine a "stability triangle" in the 
a-6 plane, for which all internal points and all points on 
the base (6=0) in the interval 0 < a < 2 correspond to a 
stable system. This triangle is shown in Fig. 5. 

The conditions for underdamped, critically damped, and 
overdamped transient response are found by inspecting the 
sign of the discrim.inant of Eq. (2.14). The resulting condi- 
tions are: 



2 

(a+6) < 46 V- - ■ > under damped (2.18) 

2 

6 1, (a+6) = 46 V- > critically damped (2.19) 

2 

a < 1, 6 < Ir (a+6) > 46 <-■ ■ > overdamped (2.20) 
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BETA (3) 




FIG. 5. THE STABILITY TRIANGLE 
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All Other values of 
(a, 6) inside the stability 
triangle and on the base 
(6=0) in the interval 
0 < a < 2 



< — V demped oscillatory natural 



mode with a rate of 
oscillation equal to one 
half the sampling frequency 



The transient response 
contains at least one 



The permissible values of a and 6 are shown in Fig. 6 



A 

2 . 




1 



6 



1 



2 



c> 



a 

Fig. 6. Allowable Values of a and 6 in a- 6 Tracker 

Further restrictions will be placed on this region, when the 
frequency response characteristics of the system are examined. 

2 . Frequency Response Characteristics 

In this section, the frequency response characteristics 
of the a -6 filter are formed. This approach has not been 
done before, in detail. 

The frequency response of the filter can be found by 

*1 (jl} * T (i\ * 

placing Z = e-* and e-' = cos wT + j sin wT into Eqs . (2.10)- 

(2.13). These equations are now in the form: 

„ (ejwT) _ [g (cos2a)T-oosgjT) +6cosajT] + ~j [a (sin2a)T-sincjT) +Ssina)T] ^2 2 
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( 2 . 22 ) 



. „ [^(cos2(jjT-cx)StoT) ] + j [;^(sin2u)T-sinu)T) ] 

= -5^ — J 



D(e^“'^) 






_ [(a+B)oosa)T - a] + j [ (g+B) sinu)T] 



(2.23) 






[|•(cos(J^)T-l) ] + j [|sino3T] 
D(e>^) 



(2.24) 



where 



D(e^'^'^) = [cos2o)T- (2-a-B) COSO)T+ (1-a) ] + j [sin2u)T- ( 2-a- B ) sincoT ] 



The amplitude and phase characteristics of G^g (smoothed 
position) and G..„ (predicted position) are plotted and shown 
in Figs. 7-16, for several values of a and B. Also the ampli- 
tude and phase characteristics of G^g (smoothed velocity) and 
G.^p (predicted velocity) are shown in Figs. 17-26. All the 
amplitudes and phases are plotted as a function of a, B and 
o)T . 

Observing these figures and Eqs . (2.7) - (2.8) one 

N N 

finds that Xg is the result of passing through a low pass 
filter, is the result of differentiating X^, and a should 
never be larger than one. 

It may be seen also, for positional smoothing, if 
a = 0, all sensor information is ignored whereas if a = 1, 
there is no smoothing of positional information. Similarly 
6=0 causes sensor information to be ignored in the estima- 
tion of velocity whereas 8 > 1 will cause overcorrection. 
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PHASE (radians) AMPLITUDE |G 
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PHASE (radians) AMPLITUDE | I 

XP XP 

-?.S0 -a. 00 -I.SO -1.00 -0.50 O.OO cP,®° 1^0 g. OO 2, .so 3.00 3.50 




u)T 




IP 

m 

'o. 00 0. 50 ^ r. 00 l‘. SO 2.00 2. SO 3. 00 ?. SO 

0)T 

FIG. 8. AMPLITUDE AND PHASE OF PREDICTED POSITION G 
FOR a = 0.1 
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FOR a = 0.5 
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PHASE (radians) AI4PLITUDE jG 





FIG. 10. AMPLITUDE AND PHASE OF PREDICTED POSITION G^p 
FOR a = 0.5 
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PHASE (radians) AMPLITUDE |G 





FIG. 11. AMPLITUDE AND PHASE OF SMOOTHED POSITION G 
FOR a = 0.75 
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PHASE G (radians) AMPLITUDE |G 





FIG. 12. AMPLITUDE AND PHASE OF PREDICTED POSITION G p 
FOR a = 0.75 
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PHASE (radians) AMPLITUDE |G 




O 

O 




FIG. 13. AMPLITUDE AND PHASE OF SMOOTHED POSITION G 
FOR a = 0.9 
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PHASE (radians) AMPLITUDE |G 

.00 -2. so -2.00 -I.SQ -l.OO -0.50 0.00 ^.00 0.05 0.90 0.95 
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FIG. 14. AMPLITUDE AND PHASE OF PREDICTED POSITION G^p 
FOR a = 0.9 
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o 

03 





FOR a = 1.5 
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PHASE G (radians) AMPLITUDE G 



O 




wT 



a 




FIG. 16. AMPLITUDE AND PHASE OF PREDICTED POSITION G^^^ 
FOR a = 1.5 
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amplitude I G 




FIG. 17. AMPLITUDE AND PHASE OF SMOOTHED VELOCITY G 
FOR a = 0.1 
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PHASE G^p (radians) AMPLITUDE |g 



in 

o 




FOR a = 0.1 
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PHASE G^g (radians) AMPLITUDE |g 




FIG. 19. AMPLITUDE AND PHASE OF SMOOTHED VELOCITY G_ 
FOR a = 0.5 
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PHASE (radians) AMPLITUDE |G 

V P 
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FIG. 20. AMPLITUDE AND PHASE OF PREDICTED VELOCITY G^p 
FOR a = 0.5 
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PHASE G (radians) AMPLITUDE [G 



O 




FIG. 21. AMPLITUDE AND PHASE OF SMOOTHED VELOCITY G^^ 
FOR a = 0.75 
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PHASE Gyp (radians) A14PLITUDE |G 



O 

<M 





FOR a = 0.75 
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PHASE G (radians) AMPLITUDE |G 





FIG. 23. AMPLITUDE AND PHASE OF SMOOTHED VELOCITY G^.g 
FOR a = 0.9 
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FIG. 25. AMPLITUDE AND PHASE OF SMOOTHED VELOCITY G^g 
FOR a = 1.5 
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FIG. 26. AMPLITUDE AND PHASE OF PREDICTED VELOCITY 
FOR a = 1.5 
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Considering the figures for several values of a < 1 , 
a seems to control the bandwidth of the low-pass filter and 
B has more control over the damping. In fact B should be 
som.ewhat smaller than a such that resonant spikes do not 
occur. 

Better accuracy occurs between smoothed and predicted 
positions, when a has very large values compared with B. 

Table I, which follows, summarizes the accuracy and 
provides an important design tool. In general the frequency 
of the input signal, the sampling time T, and the filter 
parameters a and B control the filter's response. 
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A. ACCURACY (0-1) BETWEEN SMOOTHED AND PREDICTED 
POSITION FOR WT = 0 to WT = TT 



\ ^ 
3 \ 


0.1 


0.5 


0.75 


0.9 


0.01 


1.0-0.90 


1.0-0.98 




1.0-0.99 


0.05 


1.0-0.60 


1.0-0.90 




1.0-0.94 


0.10 


1.0-0.33 


1.0-0.82 




1.0-0.89 


0.25 






1.0-0.71 




0.50 






1.0-0.50 




0.75 


' 




1.0-0.33 





B. ACCURACY (0-1) BETWEEN SMOOTHED AND PREDICTED 
VELOCITY FOR wT = 0 to 03T = TT 




0.75i 1.0 

i 

*As we expect 100% accuracy, since the vehicle was 
generally assumed to follow straight line constant 
velocity, when an a-g filter is used. 



TABLE I. ACCURACY BETWEEN SMOOTHED AND PREDICTED, 
POSITION AND VELOCITY 
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III. NOISE CHARACTERISTICS OF ct-6 FILTER-ERRORS-CRITERIA 



In the study of any filter it is essential to know the 
characteristics of the desired signals and the noise which 
excites the filter. It is desirable also to know how the 
choice of a and B affects the degree to which the noise is 
smoothed or exaggerated by the system. The description of the 
noise processes, prediction errors and methods (criteria) 
for designing an a-S filter, proceeds as follows: 

A. NORMALIZED NOISE POWER OF PREDICTED POSITION 
By making the following change of variable: 



AX = Vg T 



(3.1) 



Equations (2.7) and (2.8) describing the filter can be 
rewritten as: 



^S 


N 


- 

(1-a) (1-a) 






, N-1 


a 






= 




1 




+ 




(3.2) 


AX 




1 

-6 (1-6) i 


1 


AX 


, j 


S 

j 


1 




= [1 


1] 


'^s‘ 




(3.3) 



AX 



where now 
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(1-a) 


1 

P 

1 






a 




1 


B' 


= 




L 


(1-6) ' 






T7D 

1 



[1 1 ] 



X = 



X 



S ' 



AX 



The mean and covariance equations are defined as 






= A 



— N-1 

'X + B'x: 



(3.4) 



.N 



A ' P A ' + B • a B ' 



(3.5) 



where 



[X]^ = E[X^] , the expected value, 



P = 



P P 

^11 ^12 



P P 

21 22 



Pll = E[(X3-Xg) (X3-X3)] 



P12 = E[ (Xg - Xg) (AX -AX) ] = 



P22 = E [ (AX - AX) (AX - AX) ] 



0^.. = standard deviation of measurement error. 

XM 
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The covariance equation for the filter becomes 



^11 


N 


'■ 2 
(1-a)^ 


2(l-a)^ 


(1-a)^ 




^12 


= 


-g(l-a) 


(l-2g) (1-a) 


(1-g) (1-a) 





' ^22 



-26(l-g) (1-g) 



, N-1 



11 



12 



22 



a g 

6^ 



r 2 , 



(3.6) 



The variance of the predicted position is easily computed 
in terms of the variances from Eq. (3.6) 



2 



P + 2 P + P 
11 ^ ^12 22 



(3.7) 



The steady state solution of Eq. (3.6) is computed by leting 
P(N) equal P(N-l) and the resulting algebraic equations are 
solved; 



^11 _ 2g - 3ag + 2g^ 

2 a(4-2a-g) 

'^XM 



^12 ^ g (2a - g) 

2 a(4-2a-g) 

'^XM 
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(3.10) 



^22 _ 6(2g^ - + 2B ~ ctB) 

2 g (4 - 2g - 6 ) 

°XM 

The predicted noise power is normalized 

2 

R = (3.11) 

‘^XM 

Substituting Eqs. (3.7) - (3.10) into Eq. (3.11) yields 

a^(2 + 2S -a6) + 6(2 + a - aB) / 

a(4 - 2a + B) 

Hence, the predicted noise pov/er is plotted as a function of 
g and 6 in Fig. 27. This figure shows that by appropriately 
adjusting a and 6 the input noise power can be reduced. One 
also finds that the parameter B greatly affects the predicted 
noise power, as expected, since differentiated noise is quite 
"noisy" and B affects this quantity. 

An interesting phenomenon occurs for small values of 
g and large values of B. The noise power increases sharply. 
The reason for this may most easily be seen by looking at the 
frequency responses in Figs. 7-16. The large resonant peaks 
in the response allow a lot of noise to come through. 

B. PREDICTION ERROR DUE TO MEASUREMENT ADDITIVE NOISE ERRORS 
In the g -B tracker the observations are in the form 




X 



N 

M 



vN ^ 

X + V 
M 



N 



(3.13) 
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NORMALIZED NOISE POWER OF PREDICTED POSITION 




FIG. 27. PREDICTED POSITION NOISE POWER AS A FUNCTION 
OF a AND 3 
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where : 



N 

V 



2 

a 



V 



N M 

V V 



additive measurement noise 
for N = M 
0 for N ^ M 



and by letting 



Variance (X^ ) 



VAR(Xp'^^) 



Variance Reduction Factor = VRF 



hence 



VRF (X 



N+1 
P ' 



VAR(X 



N+1 

P 

2 



a 

V 



) 



(3.14) 



or 



VRF(X 



N+1 
P ' 



VAR(X^'^^) 




2a^+26+a3 
a ( 4-2a~3 ) 



(Steady (3.15) 
state) 



The VRF is plotted in Fig. 28 as a function of a and g. Com- 
paring this figure with Fig. 27 (normalized noise power of 
predicted position) , the two figures look alike, but this is 
not true, when both graphs are plotted together in Fig. 29. 



C. PREDICTION ERROR DUE TO CONSTANT ACCELERATION 

Constant acceleration X(t) results in fixed error in 

N+1 

predicted position X^ on reaching steady state, i.e., for 

N -V 00 . 
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NORMALIZED PREDICTION ERROR DUE TO ADDITIVE NOISE 
(VARIANCE REDUCTION FACTOR) 




FIG. 28. PREDICTION ERROR DUE TO ADDITIVE NOISE AS A 
FUNCTION OF a AND 6 
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I 

1 




NORMALIZED R, VRF 




FIG. 29. PREDICTED POSITION NOISE (R) AND PREDICTION 
ERROR DUE TO ADDITIVE NOISE (VRF) AS A 
FUNCTION OF a AND B 
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Specifically for noiseless constant acceleration trajec- 



tories 





( 3 . 16 ) 



the steady state prediction error becomes 




" 2 
XT 

6 



( 3 . 17 ) 



Common names for steady state prediction error e due to 
constant acceleration are; 

(i) Dynamic error 

(ii) Truncation error 

(iii) Systematic error 

(iv) Bias error 

Typical design procedure for balanced dynamic and random 
errors is 



D. TRANSIENT ERROR 

The transient error of the system becomes significant when 
the input is switched from one time-shared target to another 



^.Standard deviation 
^of random error 



Dynamic 

error 



that is 




( 3 . 18 ) 
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or when the target makes a sharp maneuver. Since a- 6 trackers 
are initialized by original estimates of position and velocity 
and in many applications these estimates are so inaccurate that 
the filter exhibits transient errors before settling to a 
steady-state solution. In the transient phase, when the esti- 
mates are poor, one would like to rely on the measurements 
and thus weight the errors more heavily. After several measure- 
ments have been made, the accuracy in prediction increases. 

So, the transient error yields 



X 



N+1 



I 

n=0 



(X 



N+1 



- 



T (2 - g) 



gg (4 - 2g - 6) 



(3.19) 



The transient error Dj.,, is plotted and shown in Figs. 30-32, 
as a function of g, 6 and T. Observing these figures, one 
finds, that the sampling time T greatly affects the transient 
error. For large values of T the transient error increases. 
An interesting phenomenon occurs for small values of a and 
small values of g. The transient error increases sharply. 

E. OPTIMAL STEADY STATE PLILATIONS FOR THE a-g TRACKER 
The simplest possible system is one in which g and g 
are fixed constants. The theory of such systems has been 
studied in [8]. For such systems, T.R. Benedict and G.W. 
Bordner [1] applied the calculus of variations to derive the 
following relationship: 
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TRANSIENT ERROR 



o 

n 



o 

o 



cn 

rr 

m 



C] 



+ 

S2 C^cD 
X oo 

Q 




ALPHA, a 



FIG. 30. 



TRANSIENT ERROR FOR T = 0 . 1 AS A FUNCTION 
OF a AND 3 



1 .05 
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TRANSIENT ERROR 




FIG. 31. TRANSIENT ERROR FOR T = 0.3848 AS A FUNCTION 
OF a AND 8 
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t 



TRANSIENT ERROR 




ALPHA, a 

FIG. 32. TRANSIENT ERROR FOR T = 0 . 5 AS A FUNCTION 
OF a AND 6 



r- O 



6 



(3.20) 



2 
a 

( 2 ^ 

A few years later S.R. Neal [7] used the results of linear 
estimation theory to derive the relation; 

(2a + = 83 (3.21) 

As a matter of interest, the above Eqs. are plotted and 
compared in Fig. 33. Observing these figures it is inter- 
esting to note that the relations (3.20) and (3.21) are quite 
similar for 0 ^ a ^ 0.4 and 0 ^ 3 0.1. 

Constant parameter systems suffer from the incompatible 
demands that good smoothing requires heavy damping (i.e., small 
values of a and 3 - small noise response. Figs. 27-28), while 
good response to maneuvers requires light damping (i.e., large 
values of a and 3 - small transient error. Figs. 30-32). Light 
damping and therefore large noise response, can lead to low 
probabilities of weapon sensor acquisition and problems of 
plot-to-track association. Heavy damping, implying poor 
maneuver response (large transient error) , can cause sudden 
loss of tracks (sometimes termed track death) through failure 
to associate with subsequent plots. 

These limitations led some workers to optimize for varia- 
ble parameter systems, where a and 3 are varied according to 
the state of the track. Some systems have been developed 
wherein a and 3 were initially selected arbitrarily and changed 
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.00 O.l'l 0.29 0.43 0.57 0.71 0.86 1.00 
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during program development by trial and error. Various 
operational sets of values being derived for various states 
of track. Such methods are adaptive and are usually economi- 
cal in computer use, both in terms of required storage and 
run time, but generally have no theoretically optimum adaptation. 

F. DISCOUNTED LEAST-SQUARES CRITERION 
To find linear trajectory 



Xp = X° + N T V° 



(3.22) 



. . . . N 

which minimizes sum of weighted errors between Xp and 



M 



,N-1 



.N-2 



' ' • • • ' 1 • 6 . , 



e 



N 



I 



r= 



0 



(X 



N-r 

M 



Xp)^ 



(3.23) 



where 0^6 < 1, still too many degrees of freedom for 
selecting gain terms, a and $• Minimize 






(xJJ ^-Xp"^)^6^ + 



M_2 N-2 2 

(X^ ^-X^ ^)0^ + 



(3.24) 



by solving the above equation in order to weight the differ- 
ence between measured and predicted values. This yields 
in simple gain terms 



a 



(1 - e^) 



(3.25) 



f 



p 



I 



i 







e 



(1 - 9) 



(3.26) 
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Common names for discounted least-squares a-6 filter are 
(i) Critically damped a-6 filter 
(ii) Fading-memory polynomial filter of degree 1 
Equations (3. 25) - (3. 26) are plotted as a function of (theta) 
in Fig. 34. 

More recently, processes have been developed in which 
and are made to change with time in order to continually 
compute the least-squares line through the observations. Such 
approaches assumed that errors are equally distributed in 
X and y and had a constant standard deviation. The formulae 
for changing a and 6 in this manner were worked out in [9] . 

For the incorporation of the nth measurement: 



(2M - 1) 
N (N + 1) 






6 

N(N + 1) 



(3.28) 



The above equations are plotted as a function of the number 
of measurements N, in Fig. 35. 

It is clear that, for large N, a and 6 tend to zero, 
i.e., observations will be increasingly ignored. This sug- 
gests that there should be some maximum value of N. The 
maximum value used is generally 7 to 15. Such a method may 
be made adaptive if a means of detecting changes in motion 
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FIG. 34. a- 6 GAIN FOR LEAST-SQUARE CRITERION AS A 
FUNCTION OF THETA 
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FIG. 35. a-S gain for CONTINUALLY COMPUTE THE LEAST-SQUARES 
LINE AS A FUNCTION OF OBSERVATIONS (N) 
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■L 






m 



is used, i . e. , 



if turn detection is provided. Then, if a 



turn is detected, the values of a and 6 may be raised simply 
by lowering N. Doing this will improve the turn following 
capability . 

1 . Numerical Example 

Assume the following; 



Range measurement 



X 

max 

Desire 



/vARCXp"^^) 



Op = 30.5 ft 

5g = 150 ft/sec 



31.6 ft 



Find suitable critically damped a-6 filter design. 
By substituting Eqs . (3.25) and (3.26) 

a = 1 - 



6 = (1 - e)^ 



in Eq. (3.15) 



VRF 



VAR(Xp'^^) 



and letting 



2a^ + as + 2 b 

a (4 - 2a - B ) 



VAR(X 



N+1 
P ^ 



(31.6)^, = (30.5)^, 



75 



yields 



VRF = 1.074 

0 = 0.5 

hence 

a = 0.75 

B = 0.25 

From Eq. (3.18) 

e* = 3 yVARCXp"^^) = 3 J{31.6) ^ = 94.8 ft 

but from Eq. (3.17) 

* X 

" = ■ 

hence 

T = 0.3848 sec. 



Finally, from Eq. (3.19) 




T^(2 - g) 
aB (4 - 2a - B) 
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yields 



D = D = 0.439 ft^ 

yN+1 

The values for a = 0.75 and 6 = 0.25 satisfy the critically 
damped Eq. (2.19) 

(a + 6)^ = 4 6 and 6^1 

The above values for a-6/ VRF and D can also be determined 
by using Figs. 27-31-34. 
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IV. ESTIMATING OPTIMAL TRACKING FILTER 
FOR MANEUVERING TARGET 



The work which follows is an attempt to apply optimum 
filter theory to the tactical radar environment. There it is 
desired to obtain optimal estimates of target positions and 
their predicted tracks. 

There are several possibilities for structuring a digital 
filter. From the work of Kalman [3], the filter should con- 
tain the same dynamic model as that of the incoming signal 
shown in Fig. 36. 



Measuranent Noise Vj,^ 



Excitation ^ 


Signal Process 


Signal 


Measuranent 


Observable 




Dynamics 




Process 


^K 








H 





a. Signal Process 




b. Filter 

Fig. 36. Schematic of Filter Structure in Relation to Signal Process 
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The modeling of the various system components involved 
in a tactical weapon system, such as the radar measurement 
process and the target itself, is essential in the design 
of practical tracking and control algorithm. By modeling 
the target to be tracked and the accuracy of the radar's 
measurements [10] , then a practical tracking procedure, con- 
sistent with the computer limitations and weapon system 
requirements, can be designed. 

A. SENSOR AND VEHICLE MODELING - DYNAMIC EQUATIONS 

The tracking systems under consideration utilize sensors 
that provide measurements of range and bearing. The selection 
is intended to reflect that this pair of measurements is most 
common, however, other output measurements as elevation and 
range rate (Doppler) are often available. 

The vehicles to be tracked can be modeled by the state 
equations : 



X(K+1) = <J>X(K) + Gu(K) 

where 

range at time K 
range rate at time K 
bearing at time K 
bearing rate at time K ^ 



X(K) = 



Y(K) ‘ 




Y(K) 




e (K) 




e (K) 





( 4 . 1 ) 



= vehicle state vector at time t = KT 
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'u. (K)' 




change in vehicle range rate between 


J. 




time K and time K+1 


u. (K) 




change in vehicle bearing rate between 


r- ^ 




time K and time K+1 



1 T 0 0 
0 10 0 
0 0 1 T 
0 0 0 1 



state transition matrix (-'^-•2) 



G = 



'O 0 ' 
1 0 
0 0 
0 1 



and 

T = sampling period. 



The tracking sensor measures target position, range and bearing 
or elevation and provides the following output equation: 



Y(K) = HX(K) + V(K) 



where 



measured rate at time K 
measured bearing at time K 

(4.3) 





10 0 0 


> 


'v, (K)' 


H = 




and V(K) = 


J. 




0 0 10 

4 


1 


V2(K) 
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I. 

I The measurement noise covariance matrix R(K ) , satisfies 
R(K) = E[V(K)v'^(K) ] = 

! 

assuming the noise V^(K) and V 2 (K) are independent. The 
selection of sensor coordinates (R,e), at this point, for 

' the covariance matrix R(K) has been made because the output 

i 

I matrix H assumes the extremely simple form shown, and the R(K) 
becomes diagonal. When Doppler measurements are available, 
this selection of sensor coordinates becomes extremely advan- 
tageous because the cartesian forms for H and R(K) becomes 
complex and time varying and often impose computational 
penalties for real-time implementation. 

Consider the following polar to cartesian transformation 
(i.e., from R,9 to X,Y) 



R(K) 

0 



0 (K) 



(4.5) 



X = R sin 6 



Y = R cos 9 



(4.5a) 



For the Kalman filter in XY coordinates, the measurement 
covariance matrix R(K) is a function of the radar-target 
geometry. Therefore, the elements of the covariance matrix 
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R(K) 



2 2 

^XX °XY 

2 2 

*^XY *^YY 



are 



a 



a 



a 



2 

XX 



2 

YY 



2 

XY 



2 2 2 2.2 
a cos 9 + R a sin 0 



2 2 2 2 2 
sin 0 + R Og cos 0 



2 2 2 

[Oj^ - R 0 g] sin 0 cos 0 



(4.5b) 



(4.5c) 



2 2 

where and o„ are the variances of the range (R) and 

K u 

bearing (9) measurement errors, respectively. It is inter- 
esting to note that in general, the coordinates after the 
transformation (X and Y) are not independent. The singular 
cases where they are independent occur for 0=0°, 90°, 

180°, ... or for R = — . These special cases are easily 

°0 

explained from the assumption that the measurement errors are 
Gaussian, with typical contours of equal probability, which 
are ellipses. 

Hence, if the axes lie in the directions of the axes of 
the frame of reference, the covariance terms R(K)^ 2 ' ^^^^21 
are zero and the independence of X,Y errors for 0=0°, 90°, 
etc . 

When R = o-/o„, these ellipses reduce to circles which 
may be considered as limiting ellipses whose major and minor 
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* 



axes lie along the axes of the reference frame. Finally, 
since the observation y(K) consists of x and y position 
measurements, hence, the observation V(K) can be assumed to 
have covariance 



R(K) 




(4.5d) 



B. STATISTICAL DESCRIPTION OF TARGET MANEUVER 

The input sequence u(K) is additive state (or maneuver) 
noise that results in the vehicle deviating from a constant 
velocity trajectory. 

Although the maneuver history and observation noise are 

T 

not independent, the covariance E[u(K)V (j)] is zero. Indeed, 
the radar cross section of a piloted vehicle changes during 
a maneuver, causing the radar observation noise to depend on 
the particular maneuver being exercised by the vehicle. 

The maneuver noise is neither Gaussian nor white. For 
example, the pilot of an aircraft moving at constant velocity 
will generally not maneuver unless threatened by either radar 
detection or attacking vehicles. His maneuver will then often 
be a turn or an increase or decrease in his forward velocity. 

A typical maneuver probability density is shown in 
Fig. 37. 
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Fig. 37. Typical Probability Density of Target Maneuver 

The quantity A is the maximum acceleration which the plane- 
pilot combination can withstand. Values of the density 
between non-maneuver (u = 0) and maximum (u = ±A) are non 
zero because: 

(i) The vehicle may not be accelerating at its 
maximum rate 

(ii) The projection of a circular maneuver on any 
dimension can give values of u from -A to A. 
Clearly, then, the maneuver density is not Gaussian. 

C. DISCRETE TIME EQUATIONS OF MOTION 

It is often desirable to whiten the maneuver noise so 
that system equations to which optimal filtering theory 
applies can be obtained. This is done, as in [4] . The 
whitening procedure for a discrete signal is analogous to 
the procedure developed by Wiener and Kolmogorov to white 
continuous signals [11]. So, u(k) may be expressed recursively 
in terms of the white noise sequence W(K) by 
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u(K+l) 



pu(K) + W(K) 



( 4 . 6 ) 



Using the above equation, the following set of system equations 
are obtained having white noise sequences W^^(K) and W^Ck) 
as their only inputs: 



where 



and 



X(K+1) = 4>X(K) + GW(K) 



( 4 . 7 ) 



Y(K) = HX(K) + V(K) 



( 4 . 8 ) 



X(K) = 



r (K) 
r (K) 

Ui(K)j 
0 (K) I 

e (K) 

U (K) 

^ J 



( 4 . 9 ) 



W(K) 



(K) 



W2(K) 



( 4 . 10 ) 
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Also 



1 T 
0 1 
0 0 
0 0 
0 0 
0 0 



0 0 

1 0 

P 0 

0 1 

0 0 

0 0 



0 0 

0 0 

0 0 

T 0 

1 1 

0 p. 



state Transition 
Matrix 



(4.11) 



G = 



0 

0 

1 

0 

0 

0 



0 

0 

0 

0 

0 



H = 



1 0 
0 0 



0 0 0 0 
0 10 0 



(4.12) 



= [Observation matrix] 



(4.13) 



and where 



Q(K) = E[W(K)w'^ (K) ] 



2 /I 2, 

^Ml^^ P ) 



2 M 2, 



Covariance of 
the measured 
noise 



(4.14) 
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2 

The values for the maneuver variances a... and the correlation 

Ml 

coefficient, p, depend on the maneuver characteristics of 
the vehicles being tracked. 

2 

Hence, the variances, of target acceleration are 

calculated using the model illustrated in Fig. 37. The 
target can accelerate at a maximum rate A (-A) and will do 
each with a probability and a probability of not 
accelerating at all, with an assumed uniform probability 
distribution of amplitude: 

1 - (2P, + P ) 

P(u) = ^ ^ (4.15) 



ef accelerating between -A and +A. 

The acceleration variable, therefore, has mean zero and 
variance 



2 



o 



Mi 



A^(l + 4P^ - P 2 ) 
3 



(4.16) 



Consequently, the variables u^^ and u^ / which are assumed 
independent, have zero means and u^ has variance 



2 

°M1 



A^T^ (1 + 4P^ - P^) 
3 



(4.17) 



while U 2 has variance 



2 

®M2 



A^T^ (1 + 4P^ - P 2 ) 



3R 



(4.18) 
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where R is the target range from the sensor and where all 
these quantities have appropriate units. 

The correlation coefficient (|p 1 ^1) can be modelled 
by 

_ E[u(K)u(K-l) ] 

P 2 

'^Mi 

The quantity, y, is essentially the inverse of the average 

maneuver duration. This correlation model is analogous to 

2 ~ 

that in [4], which has the discrete time form r(K) = e , 

where a is the inverse of the continuous maneuver time con- 
stant of the target. Hence, p equals e . When aT is small, 
p can be approximated closely by 1-aT, so that a and y become 
identical . 

The two extreme cases occur when p is unity and when p 
approaches zero. So, the first case represents the completely 
correlated case, and the second is the completely uncorrelated 
case, namely white noise. 

D. FILTER DESCRIPTION 

Three types of filters are considered as potential can- 
didate algorithms for tracking vehicles that are described 
by the model just discussed. These filters are the Kalman, 
the simplified Kalman and an a-g filter. 



1 - yT 



T < - 
— y 



T > - 
y 



(4.19) 
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1. The Kalman Filter 



The Kalman filter uses the augmented version of the 
model presented earlier in order to obtain white excitation 
(maneuver) noise. 

The method of computing the optimum estimate (the filter) 
is as follows: 



x(k|k) 



X(k|k-1) = (})X(K-1 |k-1) (4.20) 

X(K|K-1) + P(K|K-1)h’^(HP(K|k-1)h’^+R(K) 



. (y (K)-HX(KlK-l) ] 



(4.21) 



where 



P(K|k- 1) = d)P(K-l|K-l) (})'^+GQ(K-1 )g'^ 



(4.22) 



P(K|K) = P(K|K-1)-P(K|K-1)h'^(HP(k1k-1)h'^+R(K) ]"^HP (K|K-1) 

(4.23) 



In these equations 

(i) The double argument always denotes an estimate 

(ii) X(k|k) = Filtered estiamte of X(K) 

(iii) X(k|k- 1) = one-step predicted estimate of X(K) 

(iv) P(k|k) = Covariance of filtered estimate 

(v) P(k|k- 1) = Covariance of predicted estimate 

(vi) (y (K) -HX (K I K-1) ] = Error in the predicted 

observations 
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(vii) 



G(K) 



= p(k|k-i)h'^[hp(k1k-i)h'^+r(k) ] 



= Kalman Gain Matrix 



= Is a matrix of adjustment coefficients. 
The matrix, G(K), reflects the relative 
confidence one should have in the 
observed data as compared to the 
predicted estimate. 



The filter is initialized on the basis of two observations 
as follows: 



Yi(2) 

^Cyi(2)-yi(i)] 



0 



X(2|2) = 



(4.24) 



Y2(2) 

^Cy2 (2) -Y2 (1) ] 



0 



i 



The elements of the corresponding covariance matrix, P(2|2), 
are: 



90 



P{2\2) 



a 



2 

R 



a 




0 



0 

0 



0 



a^/T 0 

o 2 

2 20p 2 

2 2 
P'^MI ®M1 

^ Ml Ml 



0 

0 

0 



0 

0 

0 



0 

0 

0 



0 

0 

0 



2 

% 

ol/T 

0 



a'/T 



a 



M2''’ 

pa 




2 

M2 



0 

0 

0 

0 

"‘’M2 

2 

^M2 



(4.25) 

2 2 

The Quantities a,,,, should be calculated as discussed 

Ml M2 

previously (Eq. 4.17-4.18). 

A block diagram of the discrete Kalman filter is shown in 



Fig. 38. 




\\k \|K-1 ^|K-1^ 

\|K-1 " ^^-i|k-i 



Fig. 38. Block Diagram and Equations of the Discrete Kalman Filter 
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2 . The Simplified Kalman Filter 



By simplifying the maneuver model used in the Kalman 
filter, the state vector can be reduced from six to four 
elements, and the momber of independent components of the 
covariance matrix from ten to six. The model simplification 
is achieved by assuming (incorrectly) that the vehicle's 
change in velocity is uncorrelated between samples, i.e., 
the maneuver is white. 

The regular Kalman filter requires two augmented state 
variables in order to whiten the target maneuver. If the 
maneuver is assumed white, no augmentation need be performed 
and the simplification just discussed occurs. This simpli- 
fied Kalman filter can therefore also be referred to, as an 
unaugmented Kalman filter. 

The utility of this filter is greatest, therefore, 
when either sensitivity of tracking performance to assumed 
maneuver correlation is small, or when the target maneuver 
approaches whiteness relative to the sensor data rate. 

The equations for this filter and all quantities 
except the following have previously been defined: 



Q(K) 




Covariance of 
"assumed white" 
maneuver noise 



( 4 . 26 ) 



and the elements of the corresponding covariance matrix, 
P (2 I 2) , are: 
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oi/T 



P(2 2) = 



a^/T 



2 



0 

2 



o^/T 



a^/T 



2 

*^M2 „2 



(4.27) 



3 . Recursive Algorithm 

The problem is solved recursively by first assuming 
the problem is solved at time K-1. Specifically it is assumed 
that the best estimate X(K-1 |k- 1) at time K-1 and its error 
covariance matrix P (K-1 | K-1) are known. 

(i) Calculate the one-step prediction 

X(k|k-1) = (})X(K-1 |K-1) (4.28 

(ii) Calculate the covariance matrix for the one-step 
prediction 

P(K|k- 1) = ((.P (K-1 1 K-1) (l)"^ + GQ(K-1 )g’^ (4.29) 

(iii) Calculate the prediction observation 

y(K|K-l) = HX(kIk-I) (4.30) 
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(iv) Calculate the filter gain 



G(K) = P (k|k-1)h'^[HP (k|k-1)h'^+R{K) ] "^ (4.31) 

(v) Calculate the new smoothed estimate 
X(k1k) = X(K|K-1) + G(K) [Y(K)-Y(K|K-D ] (4.32) 



(vi) Calculate the new covariance matrix 

P(k|k) = P(k1k-1) - P(k|k-1)h'^[HP(K|K-1)h'^+R(K) ]"^HP(K|k-1) 

^ ^ 

or 



P(K|K) = [I - G(K)H] P(k1k-1) (4.33) 

In summary, starting with an estimate X(K-1|k-1) and its 
covariance matrix P(K-llK-l) after receiving a new observation 
Y (K) and calculating the six quantities in the recursive 
algorithm, a new estimate X(k1k) and its covariance matrix 
P(k|k) are obtained. 
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V. IMPLEMENTATION AND SIMULATION RESULTS 



In order to evaluate the three filter algorithms in a 

variety of tactical environments, an air vehicle type, and two 

tracking sensors were selected for analysis. 

Maneuver statistics (A, P^, u) were selected to match 

2 2 

the vehicle, and sensor statistics (a„, o^, T) were selected 

K U 

for each combination of sensor and data entry evaluated. 

One trajectory was constructed for the vehicle that consists 
of a straight track and a maneuvering track and is shown in 
Fig. 39. 

The x-direction vs time of range, velocity, acceleration, 
bearing, bearing rate and bearing acceleration, were plotted 
and shown in Figs. 40-45. 

For the scenario considered, an aircraft at R = 10 KM moves 

at 100 m/sec (200 knots) , can maneuver at a maximum accelera- 

2 

tion of 4g (A = 39.24 m./sec ), and has a probability of 
maneuvering at max 0.2 (2P^ = 0.2, P^ = 0.1), and a probability 
0.5 of not maneuvering at all (P 2 = 0.5). Assume a lazy 
maneuver that will provide correlated acceleration inputs for 
periods between 10 and 30 sec. Hence, with an average maneuver 
duration of 20 sec, by using Eq. (4.19), u = 1/20 = 0.05 = a, 
and the correlation coefficient P=l-aT=l- .005T. 

The two sensors were classified as: 
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FIG. 39. VEHICLE TRAJECTORY IN X-Y 
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FIG. 41. VELCOITY (X-AXIS) VS TIME 
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FIG. 42. ACCELERATION (X-AXIS) VS TIME 
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FIG. 43. BEARING (X-AXIS) VS TIME 
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FIG. 44. bearing RATE (X-AXIS) VS TIME 
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FIG. 45. BEARING ACCELERATION (X-AXIS) VS TIME 
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A. EXAMPLE 1 - AIR SEARCH RADAR 



The radar data rate was ten samples per second (T = 10 
sec) hence, the correlation coefficient p = 1 - 0.05(10) = 0.5, 
and the sensor processing noise (measurement noise variances) 
has been taken into account, 500 m in range (a = 500 m) and 

I\ 

17.4 mrad in bearing (o„ = 1 degree = 17.4 mrad) . The 
variances of maneuvering at max acceleration, and not maneu- 
vering at all, were calculated by using Eqs. (4.17-4.18) 

2 2 -4-2 

and found to be, o,,, = 46193 m/sec and = 1.346 x 10 sec 

Ml M2 

1 . Kalman Filter Evaluation 

This model can be used with the Kalman filter 
Eqs. (4.29-4.31-4.33) to determine a set of filter gains. 

For the model assumed, by using Eqs. (4.5, 4.11-4.14, 
4.25), the following matrices were obtained: 



<j>= 



1 10 
0 1 
0 0 
0 0 
0 0 
0 0 



0 0 
1 0 
0.5 0 

0 1 
0 0 
0 0 



0 0 
0 0 
0 0 
10 0 
1 1 
0 0.5 



G = 



0 

0 

1 

0 

0 

0 



0 

0 

0 

0 

0 

1 



H 



R 



1 0 
0 0 



0 0 0 0 
0 10 0 



25 



10 



0 

3.02 X 



10 



-4 
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Q 


= 


34645 

0 


, 

o 1 

o 
1 — 1 






25x10^ 


25x10 


3 0 


0 


0 


0 


25x10^ 


51193 


23096.5 


0 


0 


0 


0 


23096 


.5 46193 


0 


0 


0 


0 


0 


0 


3.02x10"^ 


3.02xlo"^ 


0 


0 


0 


0 


3.02x10"^ 


1.41x10"^ 


0.673x10"“^ 


0 


0 


0 


0 


0.673xl0"^ 


1.346x10”'^ 



2 . Simplified Kalman Filter Evaluation 

The simplified Kalman filter is achieved by assiiming 
(incorrectly) that the vehicle's change in velocity is 
uncorrelated between samples; i.e., the maneuver is white, and 
p = 0. Here, for the model assumed, by using Eqs. (4.2- 
4.4-4.5-4.26-4.27), the following matrices are obtained. 





1 


10 0 0 * 


' 


0 0 ' 






0 


10 0 




1 0 




= 








G = 








0 


0 1 10 




0 0 






0 


0 0 1 




_0 1 ^ 
















r 


1 




4 






10 0 0 


1 




25x10 


0 


H = 








R = 




— 4 




0 0 10 


1 




0 


3.02x10 
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46193 0 

0 1.346x10 





Pi* 

'25xlo'^ 


25x10^ 


0 


— 

0 


P(2l2) = 


25x10^ 


51193 


0 

-4 


0 

-5 




0 


0 


3.02x10 

-5 


3.02x10 . 

1 

-4 




0 


0 


3.02x10 


1.41x10 

j 


B. EXAMPLE 2 - 


■ SURFACE 


AND air 


SEARCH RADAR 





In this example, all quantities except the following, have 
previously been defined for the air search radar example: 
Sampling time T = 1 sec, hence p = 1-aT = 1-0.05(1) = 0.95, 
the sensor processing noise (measurement noise variances) 
has been taken into account, a = 20 m, a, = 0.1 degree = 1.74 

mrad and the variances of maneuvering were calculated and 

2 2 2 6 2 
found to be, = 461.93 m/sec , ~ 1.346x10 sec 

1 . Kalman Filter Evaluation 

For the model assumed the following matrices were 

obtained: 



*• 












1 


' 




1 


1 


0 


0 


0 


0 




0 


0 


0 

1 


1 


1 


0 


0 


0 




0 


0 


1 0 


0 


0.95 


0 


0 


0 


G = 


1 


0 


0 


0 


0 


1 


1 


0 




0 


0 


0 

1 


0 


0 


0 


1 


1 




0 


0 


- ° 


0 


0 


0 


0 


0.95 




0 


1 
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H = 



1 0 
0 0 



0 0 0 
0 10 



0 

0 



R 



400 0 

0 3.02x10 



Q 



45.04 

0 



0 

0.13x10 



P(2|2) 



400 400 0 

400 1261.93 487.83 

0 438.83 461.93 

0 0 0 

0 0 0 

0 0 0 



0 


0 


0 


0 


0 


0 


0 


0 


0 


3.02x10”® 


3.02x10”® 


0 


3.02x10"® 


7.386x10”® 


1.2787x10”® 


0 


1.2787x10”® 


1.346x10”® 



2 . Simplified Kalman Filter Evaluation 

For the model assumed the following matrices were 
obtained: 



'l 1 0 o' 




'o 


0 


0 10 0 


G = 


1 


0 


0 0 11 




0 


0 


0 0 0 1 




0 


1 















r 




1 


0 


0 


0 




400 0 


H = 










R = 


— 




_0 


0 


1 


0 




0 3.02x10 ; 

J 
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461.93 
. 0 


0 

1.346x10"^ 




400 


400 


0 


0 


400 1261.93 


0 

— 6 


0 


0 


0 


3.02x10 

— 6 


3.02x10" 


0 


0 


3.02x10 


7.386x10 



C. THE a-s FILTER EVALUATION 

The a- 6 filter considered in this paper, and used for 
simulation, is one of many varieties possible in this class, 
is more easily implemented than either the Kalman or Simpli- 
fied Kalman filters, and has been selected for evaluation 
since it is utilized extensively in tactical applications. 
Because it is designed to minimize the mean squared error in 
filtered position and velocity under the assumption of straight 
line target motion, it has little capability to track severely 
maneuvering vehicle. The a- 6 filter examined has no provision 
to adapt to different target types, as does the Kalman filter, 
since maneuver statistics are not taken into account. The 
Eqs . (2.1-2.3-2.37-2.38) for the a-6 filter evaluated are; 



N N N N 
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where 



2 (2N-1) 
N(N+1) ' 



N(N+1) 



and T = 1 sec or 10 sec 



in order to continually compute the least-squares line through 
the observations. 

D. COMPARISON OF FILTER ACCURACIES 

One hundred (100) Monte Carlo trials were made for each 
combination of tracking filter, tracking sensor and data entry 
procedure. 

Experimental determined filtered and predicted accuracies 
in vehicle range coordinates (range, range rate, and range 
acceleration) and bearing coordinates (bearing, bearing rate, 
and bearing acceleration) with means and variances of estima- 
tion error histories, were then calculated and plotted and 
shown in Appendix A. 

Table II shows a representative siiitimary of the result, 
in which the prediction accuracies of each filter were 
compared on a percentage basis to that of the Kalman filter. 

The entries in the table were determined by averaging the 
experimentally obtained percentage degradations in each of 
the above coordinates. 

The Simplified Kalman filter, and the Kalman filter 
generally performed within twenty percent (20%) of each other. 

The a- 3 filter performance, on the average, appears to 
be about equal to the above filters for the straight part 
of track, and about thirty to sixty percent (30-60%) worse. 
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Key; 1 = within 10-20 percent of the Kalman filter 
2 = within 30-60 percent of the Kalman filter 



Table II. Synopsis of the Accuracy Comparison 
of the Three Tracking Filters 



with the greatest degradation occurring for the maneuvering 
and accelerating part of track, because the gain vector quickly 
becomes too small to correct for the large estimation errors 
resulting from target maneuvers. 
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VI. SUMMARY AND CONCLUSIONS 



While no extensive analysis is implemented, it is con- 
sidered that a reasonable and unequivocal comparison of the 
filters can be made from the material presented. 

The analysis of the filters and the radar system simula- 
tions presented in this paper is considered overall to be a 
simplified but realistic model of a sophisticated system which 
could be implemented with current "state of the art" hardware. 
Based on the ensemble averages the Kalman filter obviously 
provides a somewhat better tracking response for the target 
track tested. 

The tracking ability of the a-B and Kalman filters appears 
to be about equal for "look alike" targets in close proximity, 
under the assumption of straight line motion. The a-B filter, 
however, provided unsatisfactory performance when the tracked 
vehicles executed maneuvers. Based also on the simulation 
results, the Simplified Kalman filter becomes attractive for 
implementation, because it provided tracking accuracies within 
ten-twenty percent of the Kalman filter. The utility of this 
filter is greatest, when either the sensitivity of tracking 
performance to assumed maneuver correlation is small, or when 
the target maneuver approaches whiteness relative to the 
sensor data rate. 

The filter implementation requirements increase in the 
following order: a-B filter. Simplified Kalman filter. 
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Kalman filter. Moreover, the "complexity factor" between the 
above filters is about two-to-one. 

Finally, in most applications, the answer to the question, 
"which filter is most accurate?", does not alone determine 
filter selection. Indeed, the following questions must all 
be answered in the filter selection process to obtain the 
"best" filter for a particular system: 

a. What are the actual accuracies of each filter? 

b. What are the relative filter accuracies? 

c. What are the tracking accuracy requirements of the 
system? 

d. How sensitive is system performance to tracking 
accuracy? 

e. What are the computer requirements of the filter? 

f. What are the computer limitations of the system? 

The list shows that filter selection involves careful balancing 
of filter accuracies, filter implementation requirements, 
and system performance goals and limitations. 
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APPENDIX A 



Simulation results, of experimental filtered and predicted 
accuracies, in vehicle range and bearing coordinates, with 
means and variances of estimation error histories, provided 
for comparison of each of the three filters. 

These results were generated using Monte Carlo simulation, 
with the following description of run: 
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FIG. 46. RANGE VS TIME FOR a- 6 FILTER 
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FIG. 47. VELOCITY VS TIME FOR a-6 FILTER 
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fig. 48. MEAN RANGE ERROR VS TIME FOR a-S FILTER 
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FIG. 50. VARIANCE OF RANGE ERROR VS TIME FOR a- 3 FILTER 
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FIG. 51. VARIANCE OF VELOCITY ERROR VS TIME FOR a-g 
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FIG. 52. RANGE VS TIME FOR SIMPLIFIED KALMAN FILTER 
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FIG. 53. VELOCITY VS TIME FOR SIMPLIFIED KALMAN FILTER 
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54. BEARING VS TIME FOR SIMPLIFIED KALMAN FILTER 
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FIG. 55. BEARING VS TIME FOR SIMPLIFIED KALMAN FILTER 
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FIG. 56. MEAN RANGE ERROR VS TIME FOR SIMPLIFIED KALMAN FILTER 
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FIG. 57. MEAN VELOCITY ERROR FOR SIMPLIFIED KALMAN 
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FIG. 58. MEAN BEARING ERROR VS TIME FOR SIMPLIFIED KALMAN FILTER 
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FIG. 59. MEAN BEARING RATE ERROR VS TIME FOR SIMPLIFIED KALMAN 

FILTER 
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60. VARIANCE OF RANGE ERROR VS TIME FOR SIMPLIFIED 

KALMAN FILTER 
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61. VARIANCE OF VELOCITY ERROR VS TIME FOR 
SIMPLIFIED KALMAN FILTER 



128 



(RADIANS) ^ 



.0 7.50 15. CC -22. 30.00 




C.O 7.50 15. CC 22. SC 30.00 



(SEC) 



X-Scale: 0.375 units 

Y-Scale: 0.272E - 10 units 



FIG. 62. VARIANCE OF BEARING ERROR VS TIME FOR 
SIMPLIFIED KALMAN FILTER 
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FIG. 63. VARIANCE OF BEARING RATE ERROR VS TIME FOR 

SIMPLIFIED KALMAN FILTER 
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FIG. 64. RANGE VS TIME FOR KALMAN FILTER 
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FIG. 65. VELOCITY VS TIME FOR KALMAN FILTER 
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FIG. 66. ACCELERATION VS TIME FOR KALMAN FILTER 
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FIG. 67. BEARING VS TIME FOR KALMAN FILTER 
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FIG. 68. BEARING RATE VS TIME FOR KALMAN FILTER 
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FIG. 69. BEARING ACCELERATION VS TIME FOR KALMAN FILTER 
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FIG. 70. MEAN RANGE ERROR VS TIME FOR KALMAN FILTER 
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FIG. 71. MEAN VELOCITY ERROR VS TIME FOR KALMAN 
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FIG. 72. MEAN ACCELERATION ERROR VS TIME FOR KALMAN 



22.62 



13.89 



5. 16 



-3.56 



12.29 



• 21.02 



29. 75 



FILTER 



139 



(RADIANS) 



0,0 7.50 15.00 22. 5C 30.00 




7.692E-03 



5.76PE-0'» 



-€. 535E-C3 



-1.365E-02 



-2.077E-02 



-2. 709E-O2 



-3.5 00 E-07 



X-Scale: 0.375 units 

Y-Scale: 0.0007 units 



FIG. 73. MEAN BEARING ERROR VS TIME FOR KALMAN FILTER 
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FIG. 74. MEAN BEARING RATE ERROR VS TIME FOR KALMAN FILTER 
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75. MEAN BEARING ACCELERATION ERROR VS TIME 
FOR KALMAN FILTER 
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FIG, 76. VARIANCE OF RANGE ERROR VS TIME FOR KALMAN 
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77. VARIANCE OF VELOCITY ERROR VS TIME FOR KALMAN FILTER 



144 



(M/SEC^) ^ 



0.0 7.50 15.00 22. 50 3C.C0 




1.6 19E-01 



1.3 19E-0 I 



l.Cloe-Ol 



7.190 e-02 



4. 19 IE-02 



1. 192E-02 



-1. e(7E-C2 



X-Scale; 0.375 units 
Y-Scale: 0.003 units 



FIG. 78. VARIANCE OF ACCELERATION ERROR VS TIME FOR KALMAN FILTER 
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FIG. 79. VARIANCE OF BEARING ERROR VS TIME FOR KALMAN FILTER 



146 



(RAD/SEC) ^ 



0.0 r.50 15. CC 22. 5C 30.00 




X-Scale: 0.375 units 

Y-Scale: 0.385E - 09 units 



FIG. 80. VARIANCE OF BEARING RATE ERROR VS TIME 

FOR KALMAN FILTER 
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81. VARIANCE OF BEARING ACCELERATION ERROR VS TIME 
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